Precise measurement of the VT-boson mass with the CDF II detector 
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The mass of the W boson, Mw, is an important pa- 
rameter of the standard model (SM) of particle physics. 
Precise measurements of Mw and of other electroweak 
observables significantly constrain the mass of the as- 
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yet unobserved Higgs boson, which is predicted by the 
electroweak symmetry-breaking mechanism in the SM. 
Previous measurements [l|-|4j yield a world average value 
of M w = 80 399 ± 23 MeV ||, i] and, in conjunction 
with other electroweak data, determine the Higgs boson 
mass to be M H = 89±26 GeV 0- If tne Hi gg s boson is 
observed, the comparison of its directly-measured mass 
with the SM prediction will be a powerful test of the 
model. An exclusion of the Higgs boson in the predicted 
mass range by direct searches would decisively point to 
new physics beyond the SM, for example radiative cor- 
rections from supersymmetric particles to Mvk[8|- 

The production of W bosons at y/s = 1.96 TeV at 
the Fermilab Tevatron pp collider is dominated by the 
annihilation process q'q — > W + X where X is initial- 
state QCD radiation. Leptonic decays of the W boson, 
W — > h>i {I = e,/i), provide high-purity samples that 
allow a precise measurement of Myy ■ 

In this Letter we report a measurement of M\y us- 
ing fits to three kinematic distributions in W — > [iv and 
W — >• ev decays. This measurement uses data corre- 
sponding to an integrated luminosity of 2.2 fb _1 of pp col- 
lisions collected by the CDF II detector between 2002 and 
2007, and supersedes an earlier result obtained in a subset 
of these data 0,|i]. The CDF II detector 4j is a general- 
purpose apparatus designed to study pp collisions at the 
Tevatron. In this analysis, charged-particle trajectories 
(tracks) are reconstructed and measured using a drift 
chamber (COT) Q immersed in a 1.4 T solenoidal mag- 
netic field. Electromagnetic (EM) and hadronic calorime- 



ters provide shower energy measurements as well as po- 
sition measurements via wire chambers embedded at the 
EM shower maximum. Surrounding the calorimeters, 
drift chambers [lO] identify muon candidates. Events 
are selected online if they have a muon (electron) with 
p T > 18 GeV (E T > 18 GeV) [|. 

Offline we select muon candidates defined by a COT 
track having pt > 30 GeV and associated with a 
minimum-ionizing energy deposition in the calorimeter 
and matching hits in the muon chambers. Cosmic rays 
are rejected with high efficiency using COT hit tim- 
ing [ljj. Electron candidates are required to have a COT 
track with pr > 18 GeV and an EM calorimeter cluster 
with Et > 30 GeV and must pass quality requirements 
on the COT track and the track-cluster matching. Addi- 
tionally, they must satisfy requirements on the following 
quantities: pseudorapidity (\r}\ < 1) @, the ratio of clus- 
ter energy to track momentum {E/p < 1.6), the ratio 
of energies detected in the hadronic and EM calorime- 
ters (i^Had/^EM < 0.1), and a x 2 -based difference be- 
tween the expected and observed transverse shower pro- 
files [J, Qjj]- We impose calorimeter fiduciality require- 
ments on electron candidates to ensure uniformity of re- 
sponse. When selecting the VF-boson candidate sam- 
ple, we suppress the Z-boson background by rejecting 
events with a second lepton. Events composing control 
samples of Z-boson candidates are required to have two 
oppositely-charged leptons satisfying the above criteria 
and an invariant mass (ro«) between 66 and 116 GeV 
and vector-summed px (Pt) l ess than 30 GeV. 

We define the hadronic recoil u = 1^-E; sm(9i)fij. where 
the sum is performed over calorimeter towers 113], with 
energy Ei, polar angle 6i, and transverse directions spec- 
ified by unit vectors hi. The sum excludes towers that 
contain energy deposition from the charged lepton(s). 
From pt conservation, the transverse momentum of the 
neutrino is inferred as p^ = —Pt ~ &, where p^ is the 
vector pt (Et) of the muon (electron). We calculate the 
PF-boson transverse mass as 



m T = V 2 {PtP t -Pt-Pt)- 



(1) 



To obtain high-purity samples of W bosons, we require 
30 < p T < 55 GeV, 30 < p T < 55 GeV, |u| < 15 GeV, 
and 60 < tot < 100 GeV. The final samples consist of 
470 126 (16 134) W -> ev (Z ->• ee) candidates and 
624 708 (59 738) W -> pv (Z -> pp) candidates. 

Measurements of M\y are extracted by performing 
binned maximum likelihood fits to the observed distri- 
butions of tot, Pti an d Pt using simulated line-shapes 
("templates") as a function of My/- A custom Monte 
Carlo simulation is used to generate templates between 
80 GeV and 81 GeV. The simulation includes a boson 
production and decay model, and a detailed model of 
detector response. The kinematics of W and Z bo- 
son pro duction and decay are modeled using the RES- 
BOS [l4| generator. Using the Z-boson data, we tune the 
non-perturbative form factor in resbos, which describes 
the boson pt spectrum at low pt (~5 GeV), and a s , 



which describes the boson pt spectrum at intermediate 
Pt (~15 GeV). The radiation of multiple final-state pho- 
tons is modeled with photos [H|. The photos QED 
model was checked with HORACE [la], which in addi- 
tion to a leading- logarithm calculation of multiple initial- 
and final-state photons, also performs an exact 0(a) 
calculation. We use the CTEQ6.6 [17J parton distribu- 
tion functions (PDFs) of the (anti)proton and verify that 
the MSTW2008 [3 PDFs give consistent results. The 
CTEQ6.6 and MSTW2008 PDFs yield similar estimates 
of the M w uncertainty. We quote the 68% confidence 
level (C.L.) uncertainty from the MSTW2008 ensemble 
of PDFs as a systematic uncertainty on M\y- 

The charged-lepton track is simulated using a detailed 
model of the passive material in the tracking volume 
and of individual position measurements in the COT. 
We use a highly granular lookup table to model ioniza- 
tion and radiative energy loss, multiple Coulomb scat- 
tering, and Compton scattering in the tracking volume. 
The simulation generates and propagates bremsstrahlung 
photons and conversion electrons to the calorimeter and 
includes Landau-Pomeranchuk-Migdal [l9| suppression 
for soft photon emission. Muon tracks from T, W, and 
Z-boson decays are used to determine the COT posi- 
tion measurement resolution (wl50 |im), which is imple- 
mented in the simulation as a function of radius. A helix 
fit (with beam constraint for promptly-produced tracks) 
is performed to simulate the reconstructed track. 

A high-purity sample of cosmic ray muons collected 
concurrently with the collider data is used to perform a 
precise alignment of the COT. The trajectory of each cos- 
mic ray muon is fitted to a single helix through the entire 
COT. This fit provides a robust reference for the inter- 
nal alignment of sense wires, including gravitational and 
electrostatic displacements, resulting in a 2-5 (xm preci- 
sion in relative wire positions. We remove the remaining 
weakly-constrained modes of COT deformation, based on 
the observed difference of (E/p) between positrons and 
electrons from VF-boson decays. 

We calibrate the tracker momentum scale using J/tp — > 
fifj, and T(1S) — > pp samples, by performing a maximum- 
likelihood fit of the data to simulated invariant mass tem- 
plates generated using the known mass values of these 
mesons [20j. The momentum scale is calibrated after 
alignment and energy loss corrections are derived from 
the J/tp sample. Nonuniformities in the tracker mag- 
netic field are corrected by measuring the dependence of 
the J/tp mass on the mean polar angle of the muons. The 
dependence of the momentum scale on the difference of 
the muon polar angles is used to calibrate the polar angle 
measurement and the residual bias in track curvature as 
a function of polar angle. A 4% correction to the ioniza- 
tion energy loss is applied to eliminate the dependence 
of the momentum scale on (1/pr) of the muons. 

After finalizing this calibration, we perform a mea- 
surement of the Z-boson mass in the dimuon channel 
(see Fig.QJ, initially blinded with an additive offset ran- 
domly selected from a flat distribution in the range [- 
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FIG. 1: The Z — > fi/j, (top) and Z — > ee (bottom) mass fits, 
showing the data (points), the best-fit simulation template 
(histogram) and the photon-pole contribution (shaded). The 
arrows indicate the fitting range. 



75,75] MeV. The unblinded result is M z = 91 180 ± 
12 sta t i 10 sys t MeV. This measurement is consistent with 
the world average of 91 188 ± 2 MeV @, H3, providing 
an incisive cross-check of the tracking simulation and the 
momentum scale. Subsequently, we include the Z — » [i[i 
mass measurement as a constraint on the momentum 
scale. The systematic uncertainties due to QED radia- 
tive corrections and magnetic field nonuniformity dom- 
inate the total uncertainty of 0.009% in the combined 
momentum scale. 

In the simulation of the electron cluster, nearby 
bremsstrahhmg photons and conversion electrons have 
their energies merged with that of the primary elec- 
tron. We use a custom implementation of GEANT4 [21j 
to model the distributions of electron and photon en- 
ergy loss in the solenoid coil and energy leakage into the 
hadronic calorimeter, as a function of Et and incident 
angle. Using the calibrated tracker momentum scale, we 
fit the E/p peak in W — > ev (Fig. [5]) and Z — > ee data 
in bins of Et to determine the electron energy scale and 
non-linearity of the calorimeter response. We fit the ra- 
diative tail of the E/p distribution to tune the amount of 
simulated material upstream of the COT by 2.6%. The 



EM calorimeter resolution is parametrized as the quadra- 
ture sum of a sampling term (12.6% /^Et/ GeV) and a 
constant term k = (0.68 ± 0.05)% applied to the cluster 
energy. A secondary constant term k 7 = (7.4 ± 1.8)% is 
applied only to the energies of bremsstrahlung photons 
and conversion electrons. We tune k on the width of the 
E/p peak in the W — > ev sample and k 7 on the width 
of the mass peak in a Z — > ee subsample where both 
electrons have E/p > 1.11. 

We use the tuned energy scale to perform an indepen- 
dent measurement of the Z-boson mass in the dielectron 
channel (see Fig. [I}, initially blinded with the same offset 
as used for the measurement in the dimuon channel. The 
unblinded result, M z = 91 230 ± 30 sta t ± 14 syst MeV, is 
consistent with the world average, providing a stringent 
cross-check of our EM calorimeter energy scale calibra- 
tion and electron simulation. Cross-checks of the Z — » ee 
mass measurements using exclusive subsamples consist- 
ing of electrons with E/p > 1.11 and E/p < 1.11 re- 
spectively, performed with both calorimetry and track- 
ing, give consistent results. The final determination of 
the electron energy scale combines the E/p-based cali- 
bration with the Mz measurement, taking the correlated 
uncertainty due to the QED radiative correction into ac- 
count. 
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FIG. 2: The distribution of E/p for the W — > ev data (points) 
and the best-fit simulation (histogram) including the small jet 
background. The arrows indicate the fitting range used for the 
electron energy calibration. 



The calorimeter towers containing lepton energy de- 
positions are excluded from the calculation of the recoil 
vector u. The underlying event energy in these towers 
is measured using the nearby towers in VF-boson data. 
The u resolution due to the underlying event (additional 
pp collisions) is modeled using data triggered on inelas- 
tic pp interactions (random bunch crossings). The p\- 
imbalance between the pff and u in Z — > U events is 
used to tune the recoil model, which also includes the 
response to the initial state QCD radiation and its reso- 
lution. Cross-checks of the recoil model show good agree- 
ment between TV-boson data and simulation. 

Kinematic distributions of background events passing 
the event selection cuts are included in the template fits 
with their estimated normalizations. Backgrounds arise 



Distribution 


IU-boson mass (MeV) 


X 2 /dof 


ror(e, v) 


80 408 ± 19 st at ± 18 ayat 


52/48 


Me) 


80 393 ± 21st«t ± 19s y st 


60/62 


Pr(e) 


80 431 ± 25 stat ± 22 syst 


71/62 


m T {[i,v) 


80 379 ± 16 8t at ± 16 ayat 


58/48 


Pt(m) 


80 348 ± 18 ata t ± 18 ayat 


54/62 


Pt(m) 


80 406 ± 22 ata t ± 20 ayat 


79/62 



TABLE I: Fit results and uncertainties for Mw- The fit win- 
dows are 65 - 90 GeV for the m T fit and 32 - 48 GeV for 
the p T and p 1 ^ fits. The \ 2 of the fit is computed using the 
expected statistical errors on the data points. 



Source 



Uncertainty (MeV) 



Lepton energy scale and resolution 

Recoil energy scale and resolution 

Lepton removal 

Backgrounds 

pr(W) model 

Parton distributions 

QED radiation 



2 

3 

5 

10 

4 



lU-boson statistics 



Total 



12 
19 



TABLE II: Uncertainties for the final combined result on My 



from jets misidentified as leptons, Z — > 11 decays with 
only one reconstrcted lepton, W — > tv — > £vvv, pion and 
kaon decays in flight (DIF), and cosmic rays. We esti- 
mate jet, DIF, and cosmic ray backgrounds from the data 
and Z —¥ ti and W — > tv backgrounds from simulation. 
Background fractions for the muon (electron) datasets 
are evaluated to be 7.35% (0.14%) from Z — >• it decays, 
0.88% (0.93%) from W -> tv decays, 0.04% (0.39%) from 
jets, 0.24% from DIF, and 0.02% from cosmic rays. 
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FIG. 3: The bit distribution for muons (top) and the p l T 
distribution for electrons (bottom). The data (points) and 
the best-fit simulation template (histogram) including back- 
grounds (shaded) are shown. The arrows indicate the fitting 
range. 



The fit results (e.g., Fig. [3]) are summarized in TableQ] 
As with the Z-boson mass measurements, the My/ fit val- 



ues were blinded during analysis by adding another un- 
known offset in the range [-75,75] MeV. The consistency 
of these results confirms that the VU-boson production, 
decay, and the hadronic recoil are well-modeled. System- 
atic uncertainties from analysis parameters are propa- 
gated to Mw by fitting events, generated with the param- 
eter values varied by their uncertainties, with the nom- 
inal templates. The statistical correlations between fits 
are evaluated with simulated experiments and are found 
to be 69% (68%) between rax and p T (p T ) fit values, and 
28% between p T and p T fit values. We perform a numeri- 
cal combination of the six individually fitted My/ values, 
including correlations, using the BLUE [22| method and 
obtain M w = 80 387 ± 19 MeV, with x 2 /dof = 6.6/5. 
The tot, p T an d Px fits in the electron (muon) channel 
contribute weights of 17.5% (35.5%), 13.8% (17.3%), and 
7.1% (8.8%), respectively. The systematic uncertainties 
for the combined result are shown in Table HTl 

In conclusion, we report a new measurement of the 
W-boson mass with the CDF II detector at the Fer- 
milab Tevatron using data corresponding to 2.2 fb _1 
of integrated luminosity. The measured value My/ = 
80 387 ± 12 stat ± 15 syst = 80 387 ± 19 MeV is more pre- 
cise than all previous measurements of My/ combined. 
The world average [| becomes M w = 80 390 ± 16 MeV. 
This result has a significant impact on the global elec- 
troweak fit p]]; the limit on the fitted mass of the SM 
Higgs boson has been reduced from Mh < 158 GeV to 
M H < 145 GeV at the 95% C.L. 
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